Analysis of Vibrations 
in Large Flexible Hybrid Systems* 



Olena V. Mul Delfim F. M 

olena@mat.ua.pt delfim@mat 



Torres 
ua.pt 



Department of Mathematics 
University of Aveiro 
3810-193 Aveiro, Portugal 



Abstract 

The mathematical model of a real flexible elastic system with dis- 
tributed and discrete parameters is considered, ft is a partial differen- 
tial equation with non-classical boundary conditions. Complexity of the 
boundary conditions results in the impossibility to find exact analyti- 
cal solutions. To address the problem, we use the asymptotical method 
of small parameter together with the numerical method of normal fun- 
damental systems of solutions. These methods allow us to investigate 
vibrations, and a technique for determination of complex eigenvalues of 
the considered boundary value problem is developed. The conditions, at 
which vibration processes of different character take place, are defined. 
Dependence of the vibration frequencies on physical parameters of the 
hybrid system is studied. We show that introduction of different feed- 
backs into the system allow one to control the frequency spectrum, in 
which excitation of vibrations is possible. 
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1 Introduction 

Nonlinear dynamical systems with distributed and discrete parameters are widespread 
in heavy, extractive and manufacturing industry, as well as in space-system en- 
gineering 0J|S]. Different complex dynamical processes are possible in such 
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hybrid systems, including vibrations, which always have negative influence on 
their functioning and sometimes even can result in system breakdown. 

Therefore, one of the important aims of investigation is to study possible 
vibrations in such nonlinear dynamical systems with distributed and discrete 
parameters, as spacecrafts of big size, as well as to find ways how it is possible 
to decrease harmful effect of vibrations on normal functioning of the systems. 

Spacecrafts of big size are considered, which contain some flexible elastic 
elements. Specifications for the spacecrafts (for example, in the control of form 
and accuracy of space antenna orientation) presume that intensive vibrations 
in such flexible constructions are not admissible. Passive damping of vibrations 
is often applied for elimination of vibrations, but strict specifications for the 
spacecrafts require other methods to control vibrations. One possibility is to 
use systems with active control of vibrations. The problem is to study the ef- 
fect of active control on the dynamical characteristics of the elastic construction 
with distributed parameters. There were several attempts to solve this prob- 
lem before, using mathematical models of the system which account discrete 
parameters only (see pQ and references therein). These results are not accurate 
enough. Following P, we take into account all the distributed parameters of 
the flexible elastic construction. The mathematical model becomes a nonlinear 
boundary value problem, which consists of a partial differential equation and 
nonlinear boundary conditions. 

For this new model it is possible to develop both asymptotical and numerical 
methods. This allows to obtain frequencies of possible vibrations 0. In this way 
we can also determine conditions of stability under active control of vibrations. 
After analysis of influence of different parameters of the considered systems it 
is possible to make a conclusion about the optimal parameters. 

The results of this paper may be applied in space-system engineering for 
design of new improved technical systems. 



In the given paper a flexible elastic system with distributed and discrete param- 
eters is analyzed. 

2.1 Basic Mathematical Model 

The basic mathematical model, which is used for investigations, is borrowed 
from p^: 
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where u(x, t) is the displacement of a point of the bar with abscissa x, S the 
bar cross-section area, p the unit volume mass, E the modulus of the material 
elasticity, (3 the coefficient of dissipation in the bar material, b the damping 
factor of the executive mechanism, c the rigidity of the centering spring, d the 
feedback coefficient, m the specimen mass, and I is the length of the elastic bar. 

It is necessary to point out that problem CQH© is a nonconservative bound- 
ary value problem. This fact introduces considerable difficulties for the com- 
putation of an exact analytical solution of the problem. The non-conservatism 
is caused by the presence of the odd derivatives on time both in the equation 
of motion Q and in the boundary conditions J5J) and ©■ This allows to as- 
sume that eigenvalues of the considered boundary value problem are complex 
numbers, what is in agreement with the results obtained in |51R)lllU). 

Problem is a nonclassical boundary value problem, and it is necessary 

to take into account real features of objects functioning. The aim is to obtain 
approximations for the eigenvalues of this system. We will tackle the problem 
both from a theoretical and a numerical point of view. 

2.2 Dimensionless Form of the Mathematical Model 

In order to account the influence of different system parameters onto possible 
vibrations, as well as to simplify cumbersome calculations, we introduce dimen- 
sionless variables u— j, x = y,r=y L , and dimensionless parameters 

Pa ba da m ES . . 

£l = T' » = es> v= es' ^-p-sv IT' (4) 

where a 2 = — . Then the dimensionless form of boundary problem may 
be written as follows: 

d 2 u d 2 u d 3 u 
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x = 0: u(0,t)=0, (6) 
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Mathematical model ©~0 describes dynamical processes in the system 
with distributed parameters, and such system has infinite number of degrees of 
freedom. Therefore it is natural to assume that stationary operating conditions 
in the considered system are not always single-frequency ones. 

3 Asymptotical Approach 

The first approach, which may be used for problem solving, is an asymptotical 
one. Although at present numerous efficient methods for the investigation of 
vibrations have been developed, and mathematically justified, the most efficient 
ones are the asymptotical methods, especially the averaging method 9 . This 
method is often used for the solution of boundary value problems appearing 
in the simulation of the behavior of real processes and in problems of optimal 
control. 

The idea of the averaging method is as follows: using some special operator, 
one replaces the system of differential equations under investigation by another 
system, the so-called averaged system. The last one should be simpler than the 
original system but at the same time it should describe the main features of the 
phenomenon under investigation. 

Today the averaging method is developed for various classes of differential 
equations with large and small parameters, although in most cases the problem 
of its justification remained unsolved till now. 

In case of the method of the small parameter, it is necessary to choose 
systems of equations of higher approximations whose solutions approximate the 
solutions of the original system of equations to within values proportional to 
integer powers of some small parameter. 

3.1 Basic Assumption 

In order to solve the boundary value problem ©"(01 it is expedient to apply 
the modifications of the asymptotical method of the small parameter. 

To do this, it is necessary to take into account that equations ©-0 contain 
five dimcnsionless parameters ©. These parameters, as a rule, are less than 
unit, although they have different orders of smallness. Dimcnsionless parameters 
£i and [i are connected, respectively, with the viscous friction in the material 
of the bar and the centering spring, whereas dimcnsionless parameter v is con- 
nected with the coefficient of the velocity feedback. Therefore, parameters E\, /j, 
and v are, for real dynamical systems, of magnitude order less than dimension- 
less parameters rj and S. This allow us to assume that dimcnsionless parameters 
characterizing dissipation of different physical nature are proportional to some 
common small parameter s, i.e. 

£l = ££l, fi = ep,, v — ev. 
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With such an assumption, mathematical model 0-Q may be transformed to 
the form: 

d 2 u d 2 u _ d 3 u 
dr 2 dx 2 1 dx 2 dr ' 

x = 0: m(0,t) = 0, (9) 



r d 4 u _sd 3 u 2 _ d 3 u 



dxdr 3 dr 3 dxdr 2 

d 2 u , _ _ . 9 2 w 9m „ 

+ f? ^ +e 0*j + ei) _ + _ = o. (io) 

We are now in conditions to apply for investigations some modifications of 
asymptotical methods of the small parameter. 

3.2 Conservative Case 

As it was mentioned above, mathematical model JHJ) (II Oil includes different fac- 
tors, which characterize energy dissipation in the system. Therefore both equa- 
tion of motion iJHJ and boundary conditions @ , l|lUI) contain the odd derivatives 
on time, which represent non-conservatism of the physical model of our hybrid 
system. 

But for the first approximation we can neglect all the dissipative factors of 
the mathematical model, and consider its conservative case, i.e. put e = 0. 
Then we will have E\ = ^ = v = and problem (|S|)- (|l(J[l simplifies to 

-!_--!_ = (ll) 
dr 2 dx 2 ' 1 ' 

5 = 0: u(0,r)=0, (12) 

x — 1 6 + -q + ^ — (13) 

dxdr 2 dr 2 dx 

It is easy to solve exactly problem ((TT|) - (tT3*)l and obtain the next transcen- 
dental equation for determination of the eigenvalues: 

cot(a,)= VU 2 , (14) 

where ui is an eigenvalue of the conservative boundary problem, which is a real 
number. 

Formula (|14|l can be used as the initial estimate of the complex eigenvalues 
of the original nonconservative boundary problem. 
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3.3 First Asymptotical Method 

According with [7], we will look for the solution of equation (JSJ in the following 
form: 

u{x,t) = (Acos(Ax)+Ssin(Ax))e(-^ eflA2+lA ) T , (15) 

where A is an eigenvalue of the nonconservative boundary problem, A and B 
are arbitrary constants, and i = ■ 

Solution (|15|l satisfies equation (|SJ) with accuracy e 2 . After substitution of 
solution ifTS^I into boundary conditions @ and (fTUf) , and their satisfaction with 
the same accuracy e 2 , we obtain a transcendental equation of the form 

The solution of this equation at e = is known, and it determines eigenvalues 
u) of the corresponding conservative boundary problem (|1 1ft — . 

In solving transcendental equation Hlfijl . it is necessary to account smallness 
of damping forces. This allow us to assume that the correction factor for non- 
conservatism is a quantity of order e, i.e. 

A = w + eAi, (17) 

where w is an eigenvalue of the corresponding conservative problem, and eAi is 
the correction factor for the non-conservatism of the problem. 

Substituting (|17fl into transcendental equation (|16fl . and cutting off a Taylor 
series for left-hand side of this equation as follows 

cot (w + eAi) = cot (w) — (l + cot (w) 2 ^j X t e + O (e 2 ) , 
it is possible to present the next expression for Ai: 

ir/Sw 2 [(5 p, + v 8 — E\) r]w 2 — if\ 
1 7] 2 8 2 w i + (8 r\ — 2 8 + rj) ijw 2 +7^ + 1 

Thus, eigenvalues of the nonconservative problem (|S|l- (|l(J|) are approximately 
determined in the following way: 

i e -q 8 w 2 \{8 p + v 8 — e i ) r\ w 2 — v~\ 

X~w+ — — r . (18) 

j] 2 S 2 w 4 + (8 r\ - 2 8 + rj) r]w 2 + r/ + 1 

Now, it is possible to substitute the complex value found for A into the 
exponential function of solution H15(l . Analysis of the real part of the exponent 
allows to obtain condition 



e 1 rj 2 S 2 w 4 + (2rj 2 8 2 {ii + v) + E 1 r j 2 {l - 8) - 2-q 8 e t) w 2 + £ + rj) - 2-q 8 v 



i] 2 S 2 w 4 + [77(1 + 8) — 2 8] r/w 2 + 77 + 1 



<0. 



(19) 



6 



whose satisfaction leads to self-excitation of vibrations, and self-vibrations in 
case of equality. 

It is necessary also to point out that the imaginary part of the exponent 
represents the fundamental frequency of the nonconservative system. 

Parameter v, which may be adjusted by introducing different feedbacks, en- 
ters into inequality 1(19(1 . Hence, the possibility appears to control the frequency 
spectrum, on which vibrations excitation is possible, by varying dissipative char- 
acteristics of the hybrid dynamical system. 

3.4 Second Asymptotical Method 

It is known that in a conservative system (at e = 0) periodic motions appear 
with any amplitude dependent on initial conditions. But in the nonconservative 
system (at e ^ 0) they have only fixed amplitudes, which conform to equality of 
an energy present at the expense of negative resistance and energy dissipation, 
which is always present in a real dynamical system. 

Thus, a solution of nonconservative boundary problem (|8"|l- l(10(l may be 
sought for in the form 

u(x,t) = uo (x, t) + £U\ (x,t) , (20) 

where uq (x, r) is one of the solutions of the conservative problem (|11|> — (|13|) . 
It is easy to obtain the next expression for such solution: 

Mo (x, t) = Asin(wx) cos(iut) , (21) 

where A has the meaning of amplitude, and it is necessary to select it close to 
amplitude of boundary cycle of self-sustained vibration system. 

Then, after taking into account expressions 1)20(1 and 1)21(1 . as well as rela- 
tionship ((14(1 for eigenvalues u> of conservative boundary problem 1(11(1 - 1(13(1 . the 
boundary problem (|S)l- ((lU(l may be written as follows: 



ui d u\ 



— EiAw sin(wx) sin(wr) , (22) 



dr 2 dx 2 

x = 0: ui(0,r)=0, (23) 



d 3 u 1 d 2 u\ 



dxdr 2 



dr 2 



x — 1 : Aui cos(lot) [(ijSlo 2 — 1) cos(u;) + tjuj sin(cj)] = 
+ + Auo 2 sin(wr) [{e^rjSuj 2 — 1) — fi5) cos(u) + rj5uo{D + p) sin(w)] 



(24) 

Boundary problem 1(221) - H24 J) may be considered as a mathematical model of 
the conservative system with natural frequency to, on which an external force 
of resonance frequency effects. It is known 8 that a stationary solution exists 
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only on the condition of orthogonality of external force and normal mode of a 
system. This condition determines an amplitude, which does not increase with 
time. 

Therefore, for the considered boundary problem it is expedient to choose the 
solution in the form 

u\ (x,t) = U(x) sin(u;T) , 
which reduces our problem to ordinary differential equations: 

d 2 U 

— r- + uj 2 U = -EiAj 3 sin(W) , (25) 
dx z 

x = 0: U = , (26) 

x = 1 : Aujcos(ujt) \_(r]5ui 2 — 1) cos(ti>) + ?7c<jsin(w)] 

= esin(cjr) Auj 3 t]S(D + p,) sin(w) + Aui 2 (ei (t]Suj 2 — 1) — cos(w) 

-^-{nSu^-VS-ijUuA. (27) 
ax 

Searching the general solution of equation (|25|) in the form 

U(x) = (Bi + B 2 x) cos{lox) + (Ci + C2X) sin(wa;) , 

where Bi, B2, C\, C2 are arbitrary constants, and satisfying boundary condi- 
tions (|26|l and l|27|l allows to obtain the next condition: 



C 2 ([r, 2 u; 2 (l + S) + ( V Sco 2 - I) 2 -rj\u cot (cox) - ( V 6lo 2 - l) 2 ) 
[e 1 (T ] 2 uj 2 {l-S) + (f]Suj 2 -l) 2 + r 1 ) + 2{T 1 2 cu 2 S 2 ^ + fL)~r 1 ^S)}cu ;i ~ ' [ ' 

Inequality l|28|l gives a condition for self-vibration excitation in the consid- 
ered large flexible hybrid system. In the case of equality in H28JI we have the 
most interesting case of self-excited vibrations in the system. 

The presence in (|28f) of arbitrary constant C2 signifies the dependence of 
vibrations amplitude A on the initial conditions, which were not considered in 
our boundary problem. 

Thus, it was confirmed again jS] that, in order to find some possibilities of 
vibrations control in the hybrid system, one does not need to introduce the 
initial conditions, being enough to study the frequencies of vibrations. 



4 Numerical Approach 

The other approach, which may be used for problem investigations, is a numer- 
ical one. 
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In this paper, the numerical method of the normal fundamental system of 
solutions is used for determination of the eigenvalues of the boundary value 
problem JSJ (7J|. The justification of the convergence of the method of normal 
fundamental systems of solutions, together with the theory which explains its 
good accuracy when applied to vibration problems, can be found in the mono- 
graph 3 . 

4.1 Basic Method 

Taking into account that the eigenvalues are complex numbers, the solution may 
be searched in the following form: 

u (5, r) = ( Ul (x) + iu 2 (x)) e {<i+tul)T , (29) 

where q and u> are, respectively, the real and the imaginary parts of the eigen- 
values. They are real constants, which should be determined. 

After substitution of Q29|l in the boundary problem ©^(0), the real and 
imaginary parts of the equation of motion (JjjJ are given by 

(q 2 - lo 2 ) ui - 2 quu 2 - (1 + s\q) v![ + £i uju' 2 ' = , 
2qujui + (q 2 - lu 2 ) u 2 ~ Siuju'{ - (eiq + 1) u'i = . 

In order to ensure the compactness of further expressions, let us make the 
following definitions: 

Di = v(q 2 - + 77<5<7"V + M) ~ SrjSqu; 2 ^ + //) , 

D2 = —'Sri5q 2 uo(ii + /i) + rjuu (/iSuj 2 — 2q + v8lu 2 ^j , 

^3 = S(q 2 - uj 2 ) {eifi + 77) + eirjSq (q 2 - 3lu 2 ) + q(fi5 + e x ) + 1 , 

Di = -uj(ei + fid) + EiriSuj(Lu 2 - 3q 2 ) - 25quj{ri + eifi) . 

With these notations, after separation of their real and imaginary parts, bound- 
ary conditions 10 and can be written as 



x = 0: Ul =0, 
"2 = 0; 



x = 1 : D lUl + D 2 u 2 + D 3 u[ + D A u' 2 = , 
-D 2Ul + D lU2 - Diii'i + D 3 u' 2 = . 

An application of the numerical method of normal fundamental systems 
of solutions demands a presentation of the problem as a system of ordinary 
differential equations of the first order in the normal form, which satisfy some 
boundary conditions. 
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Therefore, in order to reduce l|3U[l to the necessary form, let us enter new 
functions 

7i = u i > 72 = u 2 , 73 = u[ , 74 = u' 2 . (33) 

Then, the desired system of differential equations in the normal form may 
be written as 

'i\ = 73 , 

T"k k (34) 

7 3 = Ai7i - A 2 72 , 

.74 = ^27l + ^172 , 

where coefficients K\ and K 2 are dependent on the sought-for quantities q and 

_ g 2 - lu 2 + £l q (g 2 + lu 2 ) _ (2g + £i (q 2 + m 2 )) u 

1_ e^tu 2 + {l + e iq f ' 2 ~ £l 2 u 2 + {l + e iq f 

In new functions (|33[) . the boundary conditions (|31|l and 1321) will be trans- 
formed to the next ones: 

x = : 71 = 72 = 73 = 74 = , (35) 



x = 1 : D171 + D272 + As73 + -D474 = , 

(36) 

-D271 + -D172 - -D473 + -D374 = . 

Now it is necessary to solve the boundary value problem (|34() , l|35|l , (|36(l . By 
means of any known numerical method, one can solve the Cauchy problem for 
the set of equations (|34|l four times with the following initial conditions: 

"f jtk (a,u,q) = \ 1 l [ k .y., 3,k = 1,2,3,4, 

where the first subscript of the function 7 is a solution number, and the second 
one is a function number. 

Thus, the normal fundamental system of solutions {7j,fe(S, ui, q)} . fc=1 2 3 4 
for the set of differential equations (|34Jl may be generated. With the help of 
such system of solutions, general solution 7fe(5, w, q) of set (|34(l may be written 
as follows: 

4 

lk{x,w,q) =^Cj7j )fe (x,w,q) , 
i=i 

where Cj are some coefficients, which should be determined. 

Using the property of the normal fundamental system of solutions 

C k =7*(0,w,g), fc -1,2,3,4, 

from boundary condition l|35|l it is easy to conclude that C\ = C2 = 0. 
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Now, satisfying also boundary conditions 13tj|) . we will obtain two linear 
algebraic equations in coefficients C3 and C4: 



\C 3 E 1 {l,u; > q)+C 4 E 2 (l,cj,q)=0, 
[C 3 E 3 (l,Lj,q) + C i E 4 (l,u,q)=0, 

where coefficients E s (l, uj, q) (s = 1, 2, 3, 4) are defined by 

q) = Dl73,l(l) + -0273,2(1) + £>373,3(1) + £>473,4(1) , 

£ 2 (1, w, q) = Bi74,i(l) + £>274,2(1) + £ , 374,3(1) + £474,4(1) , 
£3(1, u,q) = -£273,i(l) +£173,2(1) -£473,3(1) + £373,4(1), 
£4(1, u, q) = -£274,i(l) + £174,2(1) - £474,3(1) + £374,4(1) ■ 
Then, the relation 



(37) 



A(w,g) = 



= 



(38) 



£i(l,o;,g) £ 2 (l,o;, q) 
£3(1,0;, q) £4(1,0;, q) 

is a necessary and sufficient condition for the existence of a non-trivial solution 
of the homogeneous set of linear algebraic equations (|37|l . 

This relation is also a necessary and sufficient condition for the existence 
of a non-trivial solution of the boundary value problem 0, ©, GJ- It is an 
equation for a pair of numbers q and to, which may be determined as its roots. 
Thus, from equation (|38|l it is possible to determine complex eigenvalues q + iu> 
of the boundary value problem. 

It is necessary also to point out that the determinant A(w, q) from expression 
l|38(l is always non-negative, which is in agreement with the results obtained in 




4.2 Improved Method 

It is obvious that the accuracy of determination of the boundary problem eigen- 
values depends on the length of the integration interval for system (I34|l . There- 
fore, in order to decrease the calculating error, it is necessary to divide the 
integration interval [0, 1] into some smaller intervals by points x = Xi, where 
i = l,n — 1, assuming xq = 0, x n = 1. 

Then solutions 7fc(x, u>,q) of boundary value problem i|34|) . (|35|l . (|36(l may 
be formed on corresponding subintervals by solutions 7^'' (cc, a;, q) of systems of 
equations 





/ M 

= 73 , 


[7f] 


/ (i) 

= 74 + 


hh 


' = K^i 


[7f] 


' = K 2l [ 



(39) 



where function 7 subscript is a solution number, and its superscript is a subin- 
terval number, i = l.n. 
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The solutions 7^ \x, Lu,q) of Ipl9"j) must satisfy boundary conditions (|3"S*j) and 
((36() and, besides, some conjugation conditions, which link them at the ends of 
the subintervals. Having in mind that points Xi of a partition of the integration 
interval [0, 1] are arbitrary, we can write the following conjugation conditions: 



[% (xi,0J,q)}' = [7fc +1) (^,w,g)]', 



(i+l)/ 



where i = 1, n — 1. In this way, solutions "fk(x,uj,q) of boundary value problem 
<|33|) . ESjl ) EHJl are defined as follows: 



jk(x,uJ,q) = J^\x,u,q), 
1, n. 



(40) 



where 

So now it is possible to use the above described technique of solution for all 
the systems of equations (139(1 . Then, accounting 1)40(1 . we get a necessary and 
sufficient condition for the existence of the non-trivial solution of the considered 
problem: 



E^XmU^q) E 2 (x n ,uj,q) 
E 3 (x n ,u,q) E4(x n ,uj,q) 



0, 



(41) 



where E\, E2, £3, E4, can be expressed by recurrence formulas. 

Again, function A(a>, q) from equation ((41(1 is always non-negative. This fact 
complicates the numerical determination of the eigenvalues and, in this connec- 
tion, we have used the method of minimum direct search for determination of 
the pairs q and u>, which minimize function A(u>, q). 



5 Discussion of the Results and Conclusions 

In this paper we have used both asymptotical and numerical methods, success- 
fully solving nonconservative boundary value problems which describe a class of 
hybrid control dynamical systems appearing in applications of space-engineering 
design. 

We present now some results of numerical calculations for the boundary 
value problem ©-Q) as well as their analysis. 

First, let us consider the results of calculations by the first asymptotical 
method. According to ((141) and 1(18(1 . the first two fundamental frequencies of 
the nonconservative system are the next ones: 

l>i = 0.3534042288, u 2 = 2.904816694 . (42) 

Now, analysis of condition 1(191) of vibrations excitation shows that a bound- 
ary value of frequency always exists in the considered hybrid system, which 
separates frequencies domains of stability and non-stability. The graph of the 
dependence of such boundary frequency on dimcnsionless parameter v is shown 
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at Figure^] where the first fundamental frequency is also present (the horizontal 
line). The graph is constructed with the following dimcnsionless parameters: 

ei = 0.005, ^ = 0.008, T) = 7, (5 = 0.1, (43) 

which take place in the real hybrid systems P^. 

The intersection point in Figure^corresponds to the value v = 0.0529760481, 
which is the maximum allowed value: for all greater values of v, the natural 
frequency lo\ is in the non-stability domain that corresponds to excitation of 
vibrations in the system. 

As it can be easily seen from Figure^] at parameters (|43ll the vibrations will 
never take place on the second frequency u>2- 

It is also possible to calculate the same boundary frequency, using the second 
asymptotical method. To do this, we use the condition of equality in i|28|) . The 
results of such calculations at the same parameters ()43(1 support the results in 
Figure n] 

Now, let us also consider, at the same parameters (|43fl . the results of calcula- 
tions of the eigenvalues of the nonconservative boundary value problem 10-10 
by the numerical method. 



According with 141|) . in Figure 2(a) and Figurc|2(b)|we present, respectively, 



the graphs of the dependence of real parts q\ and q2, of the first two eigenvalues, 
on dimcnsionless parameter v. According with (|41|l . the imaginary parts u)i, ui2 
of the first two eigenvalues are practically independent on parameter v, and 
coincide with the results l|42|) . obtained by the asymptotical method. 

It is clear that the stability domain of the system corresponds to the negative 
values of real parts q\ and qi of the eigenvalues, and when they vanish, we 
have the case of self-excited vibrations. So, from Figure 2(b) we can see that 
vibrations will not be excited at the second frequency. Besides, from Figure |2(a)[ 
it is possible to calculate the value of parameter v at which q\ = 0. This value of 
v coincides, with accuracy 10~ 3 , with the value obtained from the intersection 
point in Figure ^ by asymptotical method. 

Thus, all the results of calculations by the asymptotical and numerical meth- 
ods coincide with a sufficient accuracy. 

In this way, for each definite hybrid system it is possible to define a character 
of possible vibrations. As dimensionless parameter v is directly proportional to 
the feedback coefficient d, which may be purposefully changed by application 
of different feedbacks pQ, it is possible to change purposefully eigenvalues of 
the considered nonconservative boundary problem. This allows to avoid in the 
system an excitation of undesirable vibrations, including self-excited vibrations. 
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